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A  SPECIAL  INVESTIGATION  TO  DEVELOP  A  GENERAL  METHOD  FOR 
THREE-DIMENSIONAL  PHOTOELASTIC  STRESS  ANALYSIS 
By  M.  M.  Frocht  and  R.  Guernsey,  Jr. 

SUMMARY 


The  method  of  strain  measurement  after  annealing  is  reviewed  and 
fcund  to  be  unsatisfactory  for  the  materials  ayallabfe  In  tSHou^w 

of  the  nr  'T  f  t  iS  deacrlbed  for  the  photoelastic  determinatira  ' 

arbSraS  Sf  method  ° V  general  body  subjected  to 

-l orary  roads.  The  method  has  been  applied  to  a  sphere  sublected  tn 

accuracy!  C°"IPreSSl''e  loads'  «eults  show  posslbiUtles  o?  C 


INTRODUCTION 


general*  ihTfT  “T  plTely  Photoelastic  procedures  cannot  solve  the 

nlshef lnde”end  ^  Stf: 38  ProMe"‘-  The  photoelastic  method  fur- 

Ofthe  ftlte  of  sSJsf »tq  •  P  “hereas  the  comPlete  specification 
six  unSn-  I  lL  a  point  requires  six  relations  to  determine 

six  unknown  stress  components, 

frn  *n  °fder  to  obtain  a  sixth  relation  it  has  been  suggested  that  the 
be  madeSa^Sr  ^eS^ThS 

^  ^  enalysif ofthifsuggested  methofshowfthat fts  ^  successful"^  >  * 

?oilson's  ™SUiTth"0f 1  lals  relatively  low  values  5 

approximatelfequal  w/f*  irthS^h^^thT  ^olsson'B  ratlos 

measurement  after  annealing  breaks  down  when  STs  ratio  apprSchefl/f3111 

s- °™ri; “  ‘SStST 

obtained  by  the  numerical  integration  of  one  of  the  differential 
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equations  of  equilibrium  in  Cartesian  coordinates  rather  than  by  strain 
measurement  after  annealing  which  involves  Poisson's  ratio.  It  will  be 
shown  that  this  permits  the  determination  of  all  six  stress  components 
at  each  point  of  a  body. 

The  report  consists  of  three  parts.  The  first  part  comprises  a 
survey  and  analysis  of  the  method  in  three-dimensional  photoelasticity 
which  rests  on  the  freezing  and  slicing  processes  and  strain  measure¬ 
ment  after  annealing.  The  second  part  presents  the  theory  of  the  new 
method.  The  third  part  contains  the  application  of  the  new  method  to 
the  determination  of  stresses  in  a  diametrically  compressed  sphere. 

The  investigation  was  conducted  in  the  Photoelastic  Laboratory  of 
the  Mechanics  Department  at  the  Illinois  Institute  of  Technology  under 
the  sponsorship  and  with  the  financial  assistance  of  the  National  Advisory 
Committee  for  Aeronautics. 

The  Research  Corporation  provided  the  funds  for  the  fellowship  held 
by  Mr.  Roscoe  Guernsey,  Jr. 

Mr.  David  Landsberg,  Assistant  Research  Engineer  in  Experimental 
Stress  Analysis,  assisted  in  all  experimental  phases  of  the  work.  It 
is  a  pleasure  to  acknowledge  his  cooperation. 

Acknowledgment  is  also  due  Mrs.  Dora  L.  Frocht  for  her  assistance 
in  the  translation  of  the  paper  by  Prigorovsky  and  Preiss  (reference  1). 


SYMBOLS 


°x>  °y>  az 
Txy>  Tyz>  Tzx 

ae>  Te 


normal  stresses,  psi 

shearing  stresses,  psi 

stresses  on  an  inclined  plane,  psi 


P#  <1 

tt  ti 

p  >  q 


principal  stresses,  psi 

secondary  principal  stresses  in  XY -plane,  psi 
secondary  principal  stresses  in  XZ-plane,  psi 
normal  stresses  in  fringes 


“xy*  Hy z >  nzx 


shearing  stresses  in  fringes 
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n9y’  n0z 

e 

r,  f>  f' 

fay 


F 

F',  F" 


f 


ex>  €y>  ez 
E 


V 

P 

A 


Ac 

R 

Rc 


ax>  ay>  az 


T  .  T  .  T 

xy>  yz>  zx 


T 

ryz 


fringe  orders  at  nor*mal  incidence  for  slices  parallel 
to  XY -plane,  XZ-plane,  and  YZ-plane,  respectively 

fringe  orders  at  oblique  incidence  for  rotations  about 
Y-axis  and  Z-axis,  respectively 

angle  of  rotation  of  a  slice j  also  angle  defining  an 
inclined  plane 

isoclinic  parameters  at  normal  incidence  for  slices 
parallel  to  XY -plane,  XZ-plane,  and  YZ-plane, 
re spectively 

isoclinic  parameter  at  oblique  incidence  for  rotation 
about  Y-axis 

shear  fringe  value  of  slice,  psi  per  fringe 

shear  fringe  value  for  slices  parallel  to  XY -plane  and 
XZ-plane,  respectively 

shear  fringe  value  of  material,  psi  per  fringe  per  inch 
linear  strains,  inches  per  inch 

Young's  modulus,  psi 
Poisson's  ratio 
load,  pounds 

area  of  equatorial  plane  of  sphere,  square  inches 
area  of  surfaces  of  contact,  square  inches 
radius  of  sphere,  inches 
radius  of  contact  areas,  inches 
normal  stresses  in  terms  of  P/A 
shearing  stresses  in  terms  of  P/A 
contact  pressure  in  terms  of  P/Ac 
contact  shearing  stress  in  terms  of  P/Ac 
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SURVEY  AND  ANALYSIS  OF  EXISTING  METHODS 
Frozen  Stresses  and.  Oblique  Incidence 


Frozen  pattern.-  It  is  now  well  established  that  elastic  stress 
systems  can  be  fixed  or  frozen  into  models  made  of  certain  diphase 
plastics  and  that  such  models  with  frozen  stresses  can  be  sliced  into 
thin  sections  without  disturbing  the  fixed  pattern  (references  2,  3, 
and  4).  Observations  of  such  slices  in  a  polariscope  yield  the  rela¬ 
tive  retardations  as  well  as  the  isoclinic  parameters  at  each  point  of 
the  pattern. 

Oblique  incidence.-  The  use  of  oblique  incidence  of  a  collimated 
beam  of  light,  as  suggested  by  Drucker  and  Mindlin,  adds  materially  to 
the  information  obtainable  photoelastically  (references  5  and  6).  The 
retardation  and  isoclinics  at  normal  incidence  are  a  function  of  the 
secondary  principal  stresses  in  the  plane  of  the  slice,  while  those  at 
oblique  incidence  depend  on  the  secondary  principal  stresses  in  a  plane 
perpendicular  to  the  wave  normal  at  each  point  of  the  slice. 

The  basic  relation  for  oblique  incidence  with  rotation  about  the 
Z-axis  is  given  by  the  following  expression 

(2Fn0z)2  =  — “T-  'l[(ffx  “  az)  +  (oy  -  ax)  sin20  +  rxy  sin  2gJ 2  + 
cos2@ 

4(txz  cos  6  +  TyZ  Sin  e)2!  (1) 


The  system  of  notation  used  in  this  report  is  shown  in  figure  1. 
Normal  stresses  are  positive  when  tensile  and  negative  when  compressive. 
The  four  components  of  shear  in  the  XY -plane  are  referred  to  either  as 
the  TXy  or  the  Tyx  shear  system,  and  the  sign  of  this  system  is 
positive  when  the  shearing  components  are  as  shown  in  figure  1  (refer¬ 
ence  7,  par.  1.3).  Similarly  the  shear  system  in  the  YZ-plane  is  posi¬ 
tive  if  the  components  are  as  shown  in  figure  1.  No  signs  are  attached 
to  individual  shearing  stress  components,  their  directions  being  deter¬ 
mined  by  inspection  (reference  7>  par.  8.2). 

By  combining  the  data  from  five  stress  patterns  of  different 
obliqueness  it  is  possible  to  determine  the  three  differences  between 
the  normal  stress  components  and  the  three  systems  of  shearing  stresses 
at  each  point  in  the  slice  (for  convenience  the  plane  of  the  slice  is 
taken  as  one  of  the  coordinate  planes). 
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It  can  be  shown  that  from  the  five  quantities  obtained  with  the 
aid  of  oblique  incidence  it  is  possible  in  turn  to  obtain  the  three 
principal  shears  at  all  points  of  the  section.  This  is  equivalent  to 
determining  Mohr's  circle  for  a  three-dimensional  state  of  stress  except 
for  its  origin  which  remains  indeterminate. 

Limitations  of  purely  photoelastic  data.-  Except  for  special  cases, 
the  optical  data  by  themselves  are  insufficient  for  the  determination  of 
the  individual  principal  stresses.  This  limitation  results  from  the 
fact  that  isotropic  stress  systems  produce  no  photoelastic  effects.  Con 
sequently,  two  states  of  stress  differing  by  an  arbitrary  isotropic 
system  produce  equal  photoelastic  effects. 

The  method  employing  scattered  light,  or  the  Tyndall  effect,  which 
was  developed  in  this  country  by  Weller  (references  8  and  9)  and  inde¬ 
pendently  by  Menges  (reference  10)  suffers  from  the  same  limitation. 

The  method  of  convergent  light  employed  by  Hiltscher  (reference  11) 
and  by  Kuske  (reference  12)  makes  it  possible  to  determine  also  the 
directions  of  the  principal  stresses  but  not  their  magnitudes. 


Separation  of  Principal  Stresses 

Free  surfaces.-  The  limitation  mentioned  above  does  not  hold  at 
free  boundary  surfaces.  A  free  surface  is  subjected  to  only  two  prin¬ 
cipal  stresses,  similar  to  those  in  plane  stress  systems.  Tangential 
slices  yield  directly  the  difference  between  the  principal  surface 
stresses.  If  in  addition  a  slice  is  taken  normal  to  the  surface  and 
parallel  to  one  principal  stress,  it  is  possible  to  determine  the  indi¬ 
vidual  principal  stresses  on  the  surface  (fig.  2).  This  method  has  been 
employed  by  Leven  and  Frocht  (reference  13)  to  determine  the  principal 
stresses  on  the  surface  of  Diesel  engine  valves.  Leven  (reference  14) 
has  also  applied  this  method  to  the  problem  of  surface  stresses  in  tor¬ 
sion,  and  Hetenyi  (reference  15)  has  applied  it  to  threaded  connections. 
In  these  applications  the  faces  of  the  slice  were  oriented  to  be  normal 
to  the  direction  of  a  collimated  polarized  beam. 

In  special  cases  the  combination  of  oblique  and  normal  incidence 
leads  to  a  complete  determination  of  the  principal  stresses.  Using  this 
combination,  the  stress  distribution  in  Saint  Venant  torsion  was  deter¬ 
mined  (reference  16). 

Plane s  of  symmetry . -  For  the  special  case  where  a  plane  of  symmetry 
exists  Jessop  (reference  17)  has  developed  an  extension  of  the  Lame- 
Maxwell  equations  (Filon's  graphical  integration)  to  three-dimensional 
cases.  By  means  of  these  extended  equations  it  is  possible  to  determine 
the  stresses  along  the  axis  of  symmetry.  However,  the  method  lacks 
generality. 
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Strain  measurement  after  annealing.-  It  has  been  suggested  by  Kuske 
(reference  18)  that  mechanical  strain  measurements  after  annealing  in. 
conjunction  with  the  freezing  method  might  be  used  to  provide  the  addi¬ 
tional  relation  necessary  for  the  determination  of  the  principal  stresses 
at  a  general  point.  If  it  be  assumed  that  the  differences  between  the 
three  normal  stress  components  at  a  point  have  been  found  photoelasti- 
cally  from  equation  (l),  there  result: 


az  ~  °y  ~  czy 


a,r  -  a„  =  C 


yx 


Oz  ~  ax  -  C 


zx 


(2) 


where  the  C's  represent  constants.  If  now  the  slice  from  the  model 
containing  the  frozen  stresses  is  annealed,  the  state  of  stress  Is 
relieved  and  the  thickness  of  the  slice  at  each  point  will  return  to 
its  original  unstressed  value.  From  this  change  in  thickness,  if  it 
can  be  measured  with  accuracy,  the  strain  at  a  point  in  a  direction 
perpendicular  to  the  slice  can  be  computed.  Taking  this  direction  as 
the  Z-axis,  for  instance,  the  strain  would  be  ez.  Then  from  Hooke  s 
law*j 


ez  =  f  Ez  "  V^x  +  0yD 


(3a) 


which  may  be  written  as 


€z 


=  IE1  “  2v)Gz  +  v<az  -  °y)  +  v(°z  ~  ax)] 


(3b) 


from  which 


°z 


E6Z  -  V(jz  -  gy)  -  V(gZ  -  Ox) 
(1  -  2V) 


(4a) 


In  view  of  equations  (2)  this  gives  the  stress  component  az  after 
which  ax  o-y  are  readily  found.  The  entire  state  of  stress  has 

thus  been  determined. 
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Limitations  of  Mechanical  Strain  Measurement 

Poisson’s  ratio  equal  l/2. -  The  method  outlined  above  would  seem 
to  solve  the  problem  and  offer  a  powerful  method  of  attack.  Closer 
examination  discloses  certain  serious  difficulties.  For  most  of  the 
plastics  used  in  this  country  Poisson's  ratio  is  very  nearly  l/2  at  the 
elevated  temperatures  used  in  the  freezing  process.  If  Poisson's  ratio 
is  taken  as  l/2  then  inspection  of  equation  (4a)  shows  that  the  denomi¬ 
nator  vanishes  and  the  equation  is  quite  useless.  The  method  evidently 
fails  under  these  circumstances.  This  is  a  manifestation  of  the  general 
proposition  from  the  theory  of  elasticity  which  states  that  when  Poisson's 
ratio  is  l/2  a  given  strain  field  does  not  determine  a  unique  stress 
field,  although  the  converse  is  not  true.  A  given  stress  field  does 
determine  a  unique  strain  field. 

Thus,  when  Poisson's  ratio  is  l/2  the  general  equations  of  Hooke's 
law  may  be  written  in  the  form  below: 


crx  = 


Ov  = 


■|  Eex  +  orj^ 


3  Eey  +  a± 
2 


°z  ~  -r  ®6z  +  °i 


(5a) 


where 


CTi  = 


px  +  gy  +  oz 


For  this  set  of  equations  the  determinant  of  the  coefficients  vanishes 
and  the  matrix  is  of  rank  two.  There  are,  therefore,  infinitely  many 
solutions.  Furthermore,  these  solutions  differ  from  each  other  by  arbi¬ 
trary  isotropic  stress  systems  gi# 

As  already  stated,  in  the  method  employing  strain  measurement  after 

annealing  the  procedure  is  to  obtain  photoelastically  two  expressions 

giving  the  differences  between  the  normal  stress  components  and  to  combine 

these  with  one  expression  based  on  the  measurement  of  one  strain  €  nor- 

z 

mal  to  the  face  of  the  slice.  This  yields  the  following  equations: 


t 
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(5b)  * 


As  in  the  case  of  equations  (5a)  the  determinant  of  the  coefficients  of 
equations  (5b)  also  vanishes  and  the  matrix  is  of  rank  two.  There  axe, 
therefore,  inf initely  many  solutions  to  equations  (5b).  Consequently, 
the  method  of  strain  measurement  after  annealing  fails  when  V  =  1/2. 

The  situation  is  exactly  analogous  to  that  existing  in  photoelasticity 
'  itself  where  the  principal  stresses  at  a  point  can  be  determined  only 
up  to  an  arbitrary  isotropic  system.  Strain  measurement  when,  v  =  1/2 
thus  adds  nothing  to  the  general  solution  of  the  three-dimensional 
problem. 

Poisson's  ratio  nearly  l/2.-  If  Poisson's  ratio  is  slightly  less 
than  1/2,  the  method  of  strain  measurement  after  annealing  should  theo¬ 
retically  lead  to  a  solution.  However,  other  difficulties  arise.  Equa¬ 
tion  (4a)  may  be  written  as 


ECZ  ~  v(Czy  +  ^zx) 
(1  -  2V) 


(4b) 


§ 


If  V  is  the  true  value  of  Poisson's  ratio  and  V  +  Av  is  its  experi 
mentally  determined  value,  then  the  error  Acxz  in  the  computed  stress 
cjz '  =  crz  +  Acrz  for  a  measured  ez  is 


tfz  *  -  crz  —  Aaz  — 


Av 


1  -  2v  -  2AV 


(erx  +  ay) 


(6a) 


With  V  very  close  to  l/2,  a  very  small  error  AV  in  Poisson's  ratio 
may  lead  to  large  errors  Aaz  in  the  computed  stress.  For  instance, 
if  v  =  0.48  is  assumed,  which  is  the  approximate  value  for  Fosterite, 
and  Av  is  taken  as  only  0.01,  then 

Aaz  =  ^”(nx  +  ffy)  =  0-50(ax  +  cry) 


(6b)  • 
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Experimental  measurements  of  v . -  Experiments  with  the  determination 
of  Poisson’s  ratio  for  Fosterite  and  similar  plastics  indicate  that  it 
will  he  rather  difficult  to  determine  the  value  of  V  closer  than 
±5  percent.  The  error  Aaz  in  equation  (6b)  would  be  particularly  large 
when  ctx  and  a-y  happen  to  be  of  the  same  sign  and  each  is  large  in 
comparison  with  az.  All  things  considered,  no  great  accuracy  can  be 
expected  from  this  method  so  long  as  Poisson's  ratio  is  nearly  l/2. 

It  must  be  pointed  out,  however,  that  strain  measurements  may  serve 
a  useful  purpose.  Assuming  that,  in  some  way  or  other,  the  normal  stress 
components  have  been  found,  the  strains  can  be  calculated  and  compared 
with  those  found  experimentally.  Here  the  error  in  the  computed  strains 
due  to  an  error  in  Poisson's  ratio  is  given  by 

A*z  =  -|(ax  +  ay)Av  (7) 


which  is  not  large. 


Method  Suggested  by  Prigorovsky  and  Preiss 

The  method  outlined  above  for  the  separation  of  the  principal 
stresses  which  employs  oblique  and  normal  incidence  of  collimated  polar¬ 
ized  light  and  strain  measurement  after  annealing  is  not  the  only  possi¬ 
ble  procedure.  Prigorovsky  and  Preiss  suggest  two  alternative  methods 
in  reference  1.  Their  procedures  combine  (l)  stress  patterns  from  nor¬ 
mal  and  oblique  incidence  with  (2)  axis  patterns  from  convergent  polar¬ 
ized  light  and  (3)  strain  measurements  after  annealing.  The  significant 
point  lies  in  the  fact  that  the  final  equations  obtained  from  their  pro¬ 
cedures  have  identically  the  same  limitations  as  equations  (5). 


THEORY  OF  SHEAR  DIFFERENCE  METHOD 


General  theory.-  A  method  for  determining  stresses  in  three- 
dimensional  problems  is  now  proposed  which  is  completely  general.  With 
this  method  the  six  stress  components  at  any  point  may  be  found.  It  is 
essentially  an  extension  to  three  dimensions  of  the  method,  long  and 
effectively  used  for  plane  problems,  which  is  known  as  the  shear  differ¬ 
ence  method  (reference  7>  ch.  8). 
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Consider  an  arbitrarily  loaded  unsymmetrical  model  with  the  set  of 
coordinate  axes  as  shown  in  figure  3*  Let  a  straight  line  AB  be  drawn 
through  point  (i)  from  boundary  to,  boundary  and  let  this  line  be  taken 
as  the  X-axis.  At  any  point  along  this  line  the  first  differential  equa¬ 
tion  of  equilibrium,  with  body  forces  neglected,  is 


dgx  +  dTyx  +  dTzx 
8x  8y  8z 


(8) 


and  upon  integration  the  stress  at  any  point  (j)  is  given  by 


where  (crx)a  denotes  the  stress  at 
any  other  point  (j)  on  the  line  AB. 

rate  of  change  of  Tyx  with  respect  to  y  and  • — —  is  the  rate  of 

change  of  with  respect  to  z.  Thus  if  values  of  were  computed 

along  a  line  through  (i)  parallel  to  the  Y-axis  and  the  curve  Tyx  =  f(y) 


point  A  and  (tfx)jj  the  stress  at 

The  partial  derivative  — is  the 

8y 


were  plotted,  then 
point  (i).  Similarly 


fid 
fir) 


would  be  the  slope  of  the  curve  f(y)  at 

is  the  slope  of  the  curve  Tzx  =  g(z) 
i 


at  point  (i).  As  in  the  plane  problem,  these  slopes  may  be  approximated 
by  computing  the  shearing  stresses  at  points  near  (i)  on  opposite  sides 
of  the  line  AB  and  forming  finite  difference  quotients.  Thus,  choosing 
neighboring  points  C  and  D  in  the  XY -plane  and  similarly  points  E  and  F 
in  the  XZ -plane. 


/8Tyx' 

1  .  (Tyx)c  -  (Tyx)c  / 

^ATyx\ 

V  8yy 

1 

1 

•H 

(10a) 
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Thus  if  the  shearing  stresses  can  he  determined  along  four  auxiliary 
lines,  parallel  to  and  on  opposite  sides  of  AB,  two  lines  in  the  XY -plane 
and  two  in  the  XZ -plane ,  one  has  all  the  data  necessary  to  obtain  the 
quotients  on  the  right  side  of  equations  (10)  and  hence  good  approxima¬ 
tions  to  the  partial  derivatives.  In  evaluating  equations  (10),  care 
must  be  taken  to  attach  the  proper  signs  to  the  shear  systems  Tyx  and 
Tzx>  as  in  figure  1. 

Substituting  the  above  approximations  for  the  partial  derivatives 
in  equation  (9)  and  replacing  the  integrals  by  summations,  the  following 
equation  is  obtained: 


<3  Ar-u-y  J  . 

K)j  -  K)a  -  -  SI 

a  a 


(lla) 


The  summations  are  evaluated  graphically  in  the  same  manner  as  in  plane 
problems.  For  convenience.  Ay  and  Az  may  be  taken  equal  to  Zac. 
Then  equation  (lla)  becomes 

J  A 

(0x)j  =  (0x)a  ±  'S  ^Tyx  ±  5~Atzx  (lit1) 

a  a 


in  which  At  „  and  At  have  the  mean  value  in  each  interval  Zac. 

Shearing  stresses  in  first  slice.-  In  order  to  carry  out  this 
integration,  it  is  necessary  to  determine  the  magnitude  and  direction 
of  the  shearing  stresses  t^  along  the  two  auxiliary  lines  in  the 
XY -plane  and  of  Tzx  along  the  two  auxiliary  lines  in  the  XZ-plane. 
The  shearing  stresses  Ty*  are  obtained  from  a  slice  in  the  XY -plane 
containing  AB  in  its  middle  surface.  The  stress  pattern  of  this  slice 
from  normal,  incidence  will  give  the  difference  between  the  secondary 
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principal  stresses  in  the  plane  of  the  slice  at  all  points,  and  the 
corresponding  isoclinics  furnish  their  orientation.  The  magnitude  of 
the  shearing  stress  Tyx  at  any  point  will  then  be  given  by 


Tyx  =  |(p*  -  q’)  sin  20 1 


(12a) 


where  p'  and  q'  are  the  secondary  principal  stresses  in  the  XY -plane 
and  0'  is  the  isoclinic  parameter.  The  directions  are  determined  by 
inspection  as  in  paragraph  8.2  of  reference  7-  Using  equation  (12a), 
the  shearing  stresses  Tyx  along  the  auxiliary  lines  and  along  AB  itself 
may  be  found. 

Shearing  stresses  in  second  slice.-  A  second  slice  lying  in  the 
XZ-plane  and  containing  line  AB  would  furnish  similar  information  for 
t •  Here  a  practical  difficulty  arises  since  the  first  slice  removes 
an  essential  part  of  the  second  slice.  One  of  several  procedures  may 
be  used  to  eliminate  this  difficulty. 

(l)  In  the  general  case  two  identical  models,  identically  loaded, 
may  be  used,  one  for  the  XY  slice  and  one  for  the  XZ  slice.  The  shearing 
stresses  Tyx  for  the  XY  slice  are  calculated  from  equation  (12a).  Simi¬ 
larly,  the  shearing  stresses  T zx  for  the  XZ  slice  are  given  by 


Tzx  =  |(p"  -  q")  sin  20"  (12b) 

where  p"  and  q"  axe  the  secondary  principal  stresses  and  0"  is  the  • 
isoclinic  parameter  in  the  XZ  slice . 

(2)  In  large  models  it  may  be  possible  to  use  a  sub  slice  from  the 
main  slice  for  determining  t^.  After  the  data  are  obtained  from  the 
main  slice,  a  small  section  containing  the  line  AB  is  cut  from  it,  as 
shown  in  figure  4(a).  The  necessary  data  in  the  XZ-plane  are  then 
obtained  from  normal  incidence  on  the  sub  slice,  as  shown.  This  procedure 
is  feasible  if  the  model  is  large  so  that  the  main  slice  can  be  made  of 
sufficient  thickness. 

(3)  In  the  particular  case  where  a  plane  of  stress  symmetry  exists, 
such  as  the  XY -plane  in  figure  4(b),  advantage  can  be  taken  of  this 
symmetry.  Referring  to  figure  4(b),  let  it  be  assumed  that  the  stresses 
on  line  AB  are  required.  The  first  slice  is  made  parallel  to  the 
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XY -plane  and  contains  the  line  AB  as  shown  in  figure  4(c).  The  orthogo¬ 
nal  slice  is  cut  from  the  opposite  side  of  the  body  so  that  it  contains 
the  symmetrically  placed  line  A'B',  along  which  the  stresses  are  the 
same  as  along  AB  itself.  The  necessary  shearing  stresses  are  calculated 
as  outlined  in  procedure  (l),  above. 

Normal  stresses.-  The  starting  value  (ax)a  will  he  obtained  from 
boundary  conditions  and  boundary  fringe  orders.  The  integration  may 
then  be  carried  out  and  values  of  ax  obtained  along  AB.  Further,  from 
Mohr's  circle  or  other  considerations: 

(crx  -  ay)j  =  (p'  -  cj')j  cos  20'j  =  2F'n'j  cos  20 ’ ^  (13a) 


(ax  -  az)j  =  (p"  -  q")^  cos  2(f)" ^  =  2F"n'\  cos  2 0"j  (13b) 


where  the  F's  and  n's  .  denote,  respectively,  fringe  value  of  the  model 
in  shear  and  fringe  order  at  point  (j).  From  equations  ( 13a)  and  (13b) 


(tfy)j  =  (ax)j  -  (p '  -  q')j  cos  20' j  (13c) 

(tfzh  =  (crx).  -  (p"  -  q").  cos  20".  (13d) 

J  J  J  J 

All  the  necessary  data  for  the  evaluation  of  (tfy)j  an(^  (°z)j  are 
obtained  from  the  slices  in  the  XY-  and  XZ-planes,  respectively. 

Use  of  oblique  incidence.-  At  this  stage  five  of  the  six  stress 
components,  namely  crx,  cry,  crz,  Tyx.>  and  tzx,  have  been  found  at 
all  points  of  AB.  There  remains  one  unknown  stress  component  TyZ.  The 
shear  system  t yZ  has  no  influence  on  the  stress  patterns  from  normal 
incidence  for  either  one  of  the  slices  but  will  have  an  effect  on  the 
patterns  from  oblique  incidence . 

In  order  to  find  the  shear  system  TyZ  an  oblique  stress  pattern  is 
obtained  from  either  one  of  the  two  slices.  For  concreteness  assume  that 
the  slice  parallel  to  the  XY -plane  is  used  and  that  it  is  rotated  in  a 
clockwise  direction  about  the  Y-axis  through  an  arbitrary  angle  0. 


Figure  5  shows  a  view  of  a  small  element  as  seen  from  the  positive  end  of  the  Y-axis.  From 
equation  (l)  the  fringe  order  n@y  at  any  point  due  to  the  oblique  incidence  is  given  by 
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and 


sin  200y  = 


2(nyZ  cos  6  +  nyx  sin  0) 
n©y  cos  Q 


(15b) 


Similar  equations  may  be  written  for  rotation  of  the  slice  parallel  to 
the  XZ-plane  about  the  Z-axis. 

It  is  to  be  noted  that  in  general  the  retardation  observed  at  any 
point  depends  on  the  direction  of  rotation  of  the  slice.  For  one  direc¬ 
tion  of  rotation  the  'fringe  order  at  a  point  will  be  different  from  its 
value  for  the  other  direction.  In  the  particular  case  when  the  slice 
contains  a  principal  plane,  then  Tzx  =  Tzy  =  0  and  equations  (l4a)  and 

(l4b)  become  identical.  In  such  cases  the  direction  of  rotation  is  imma¬ 
terial.  In  dealing  with  general  slices  it  is  important  to  note  carefully 
the  direction  of  rotation  relative  to  the  wave  normal  and  to  attach  the 
proper  signs  to  all  the  stresses. 

Each  of  equations  (l4a),  (l4b),  (15a),  and  (15b)  may  be  solved  for 
the  unknown  shear  system  nyZ.  It  is  necessary  only  to  determine  the 
fringe  order  and  the  isoclinic  parameter  0@y  along  the  line  AB.  If 
the  rotation  is  counterclockwise  equation  (l4b)  or  (15b)  is  appropriate. 
Of  these,  equation  (15b)  is  much  the  simpler.  Using  equation  (15b)  and 
6  =  45°  there  is  obtained 

xiyZ  =  |  ney  sin  20ey  -  nyx  (15c) 

With  this,  nyZ  '  and  Tyz  =  2FriyZ  are  easily  computed.  All  six  com¬ 
ponents  of  stress  are  thus  determined  for  the  point  (i),  and  therefore 
the  principal  stresses  themselves  are  determined  at  the  given  point. 

Stresses  in  the  plastic  range.-  It  should  be  noted  that  the  method 
described  in  this  report  is  not  limited  to  the  elastic  state.  It  is 
equally  valid  in  the  plastic  state  of  the  model,  provided  the  stress- 
optic  law  for  this  range  is  known.  This  follows  from  the  fact  that  the 
only  equations,  in  addition  to  the  stress-optic  law,  used  in  this  method 
are  the  equations  of  equilibrium,  which  are  independent  of  stress-strain 
relations.  It  should,  however,  be  noted  that  the  results  will  apply  to 
the  model  only  and  that  they  are  not  directly  transferable  to  the 
prototype . 


1 6 


NACA  TN  2822 


Effect  of  Poisson's  ratio.-  In  conclusion,  it  should  be  noted  that 
in  transferring  the  results  from  three-dimensional  photoelastic  models 
to  metal  prototypes  the  effect  of  Poisson’s  ratio  will  have  to  be  con¬ 
sidered.  It  is  fortunate,  as  shown  by  the  theoretical  solutions  obtained 
to  date,  that  the  effect  of  Poisson’s  ratio  on  the  most  significant  . 
stresses  is  small  (references  19  and  20). 


APPLICATION  OF  SHEAR  DIFFERENCE  METHOD  TO 
A  DIAMETRICALLY  COMPRESSED  SPHERE 
Description  of  Apparatus 

The  apparatus  used  in  this  investigation  consisted  of  the  following 
items : 

(1)  An  electric  furnace  with  temperature  controls  and  built-in 

loading  frame  with  special  jigs 

(2)  An  8-inch  photoelastic  polariscope  with  a  special  immersion 

tank 

(3)  An  oblique  incidence  jig 

(4)  A  Babinet-Soleil  compensator 

A  photograph  of  the  electric  furnace  is  shown  in  figure  6.  This  is  a 
relatively  large  oven  46  inches  high,  42  inches  wide,  and  19  inches  deep. 
It  is  fitted  with  automatic  temperature  controls  by  means  of  which  any 
desired  thermal  cycle  could  be  imposed  on  the  model.  The  furnace  was 
equipped  with  a  built-in  loading  frame  suitable  for  the  application  of 
all  basic  types  of  loading. 

A  special  loading  jig  built  for  the  investigation  is  shown  in  fig¬ 
ure  7*  It  consisted  of  a  smooth  circular  shaft  about  l/2  inch  in  diam¬ 
eter  passing  through  a  pair  of  smooth,  lubricated  guide  holes  carefully 
alined  so  that  the  axis  of  the  shaft  was  perpendicular  to  the  base.  The 
load  was  applied  to  the  top  of  the  shaft  through  a  hard  steel  ball.  This 
jig  was  found  to  give  almost  perfect  vertical  loading  and  the  friction 
was  negligible. 

A  special  jig  was  also  built  for  oblique  incidence.  The  frame  of 
the  jig  could  be  rotated  about  a  vertical  axis  through  any  desired  angle 
which  could  easily  be  measured  to  one-tenth  of  a  degree.  The  slice  was 
mounted  in  the  frame  of  the  jig  and  the  whole  unit  was  placed  in  an 
immersion  tank  with  a  suitable  mixture  of  Halowax  and  mineral  oil. 
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The  remaining  equipment  was  standard  apparatus  in  photoelastic 
laboratories  the  descriptions  of*,  which  are  available  in  the  literature. 


Test  Procedure 

Model . -  The  sphere  was  machined  from  a  cylinder  of  Fosterite  which 
was  previously  annealed  to  remove  initial  stresses.  The  machining  was 
carried  out  in  a  lathe.  The  rough  cylinder  of  Fosterite  was  first  turned 
to  a  true  cylinder.  In  order  to  form  the  sphere  a  tool  bit  was  set  in 
a  special  jig  riding  on  the  carriage  of  the  lathe.  This  bit  could  be 
rotated  about  a  vertical  axis  lying  in  the  plane  of  the  lathe  centers. 

The  cutting  was  performed  by  swinging  the  tool  bit  by  hand  around  its 
vertical  axis  while  the  cylinder  was  rotating,  and  the  radius  of  the  * 
sphere  was  slowly  reduced  by  bringing  the  tool  bit  gradually  closer  to 
its  axis.  In  this  way  it  was  possible  to  shape  the  complete  sphere 
except  for  a  relatively  small  nipple  near  the  chuck.  The  final  diameter 
of  the  sphere  was  3.313  ±  0.002  inches. 

Loading. -  The  sphere  was  placed  in  the  loading  jig  and  carefully 
alined  for  diametral  compressive  loading.  The  model  was  then  heated  to 
162  F  in  the  electric  furnace,  the  rate  of  heating  being  about  10°  F 
per  hour.  A  load  of  172  pounds  was  then  applied  to  the  model.  After  a 
soaking  period  of  about  2  hours  the  temperature  was  lowered  at  the  rate 
of  4°  F  per  hour  to  room  temperature.  The  final  diameter  of  the  equa¬ 
torial  plane  was  found  to  be  3*33^-  inches  and  the  load  axis  was  measured 
as  3.102  inches.  Although  relatively  large  deformations  were  developed 
in  the  loaded  regions,  the  main  body  of  the  sphere  was  not  notably  dis¬ 
torted  from  its  original  shape.  The  stress  pattern  of  the  whole  sphere 

in  figure  8  shows  that  the  loads  and  the  stresses  were  rotationally 
symmetrical. 

Slicing.-  In  preparation  for  slicing  the  center  lines  of  all  the 
slices  were  carefully  scribed  on  the  sphere  using  the  flat  spots  in  the 

loaded  regions  as  datum  planes.  The  slices  were  then  sawed  out  roughly 

on  a  band  saw  to  a  thickness  of  about  3/8  inch.  They  were  subsequently 
ground  by  hand  to  about  l/8-inch  thickness  in  most  cases.  Great  care 
was  taken  to  keep  the  slices  symmetrical  with  respect  to  their  center 
lines. 

Figure  9  shows  the  slicing  plan.  The  first  slice  removed  was 
parallel  to  the  equator  and  midway  between  the  equator  and  the  load 
point.  Then  from  the  opposite  side  of  the  sphere  a  meridian  slice  was 
removed.  Next  a  slice  containing  the  equatorial  plane  was  cut.  Finally 
a  slice  parallel  to  the  meridian  slice  and  halfway  out  on  the  radius  was 
removed. 
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Stress  patterns  and  isoclinics.-  The  slices  were  mounted  in  the 
oblique  incidence  jig  and  stress  patterns  at  normal  and  oblique  incidence 
were  recorded  photographically.  Typical  stress  patterns  are  shown  in 
figures  10  to  15.  Most  of  the  normal  incidence  patterns  show  very  few 
fringes.  In  order  to  obtain  accurate  data  in  these  cases  a  Babinet- 
Soleil  compensator  was  used  to  obtain  the  fringe-order  distribution 
along  the  lines  of  interest  by  point -by-point  exploration.  A  black 
cross  was  made  on  the  screen  to  mark  the  point  on  which  attention  should 
be  centered.  The  model  slice  was  then  adjusted  until  the  line  of 
interest  on  the  image  ran  true  on  the  intersection  of  the  cross  as  the 
straining  frame  was  moved  laterally.  Then  beginning  at  the  outer  edge 
and  moving  the  straining  frame  by  a  known  amount  after  each  observation 
the  fringe  order  was  obtained  at  a  series  of  points  along  the  line. 

From  these  data  the  curve  of  fringe-order  distribution  could  be  plotted. 

.  It  may  be  noted  that  fractional  fringe  orders  can  also  be  obtained 
by  the  Tardy  method  of  compensation,  the  accuracy  being  comparable  wi 
that  of  the  Babinet-Soleil  compensator. 

Isoclinic  lines  were  recorded  by  one  of  two  methods.  For  the  over¬ 
all  picture  the  isoclinic  lines  were  photographed  in  most  cases  (figs,  lb 
and  17).  From  the  photographs  averaged  sketches  were  prepared  and  used 
in  making  the  calculations.  On  several  lines  direct  sketching  of  the 
isoclinics  was  used  with  attention  being  confined  to  the  particular  line 
of  interest.  Here  the  intersections  of  successive  isoclinic  lines  with 
the  line  of  interest  were  obtained  at  very  short  intervals,  from  which 
the  distribution  of  the  isoclinic  parameters  along  the  line  could  e 
plotted.  This  method  was  found  to  be  accurate  and  considerably  less 
time  consuming  than  the  photographic  method.  White  light  was  used  m 
all  isdclinic  work. 

In  plane  stress  systems  the  isoclinic  parameter  at  a  point  on  a 
free  boundary  is  determined  by  the  tangent  to  the  boundary  at  the  poin  . 
The  isoclinic  parameter  thus  changes  from  point  to  point  along  the 
boundary  in  general.  This  is  not  necessarily  true  for  isoclmics  of 
secondary  principal  stresses.  In  the  slice  parallel  to  the  men  lan  e 
secondary  principal  stresses  at  the  boundary  consist  solely  of  one  normal 
stress  crz  which  is  horizontal.  The  boundary  is  therefore  a  zero 

isoclinic  and  no  other  isoclinic  may  intersect  the  boundary  at  any  point. 
The  higher-order  isoclinics  therefore  all  lie  within  the  boundary  forming 
closed  loops  in  this  case,  as  shown  in  figure  17- 

Fringe  value. -  The  fringe  value  of  the  material  was  obtained  from  a 

small  cylinder  about  l/2  inch  in  diameter  and  l|  inches  long.  This  was 

loaded  in  compression  in  the  special  loading  jig  used  for  the  sphere  and 
subjected  to  the  same  thermal  cycle  as  the,  sphere.  A  portion  of  the 
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cylinder  was  machined  away  to  leave  a  V  shape  as  shown  in  figure  18. 

The  resulting  stress  pattern  was  then  photographed  (fig.  19).  The 
V  shape  was  used  to  make  clearly  visible  the  fringe  of  zero  order 
occurring  at  the  sharp  edge  of  the  wedge.  In  the  cylinder  itself  the 
first  few  fringes  crowd  together  near  the  boundary  of  the  cylinder  and 
it  is  practically  impossible  to  identify  the  zero  fringe.  From  the 
stress  pattern  in'  figure  19  it  was  a  simple  matter  to  plot  fringe  order 
against  position,  which  for  the  wedge  described  is  a  straight  line 
(fig.  20). 

In  this  way  the  fringe  order  at  the  point  where  the  plane  faces  of 
the  V  intersect  the  curved  boundary  of  the  cylinder  could  be  accurately 
determined.  At  the  same  point  the  thickness  h  could  also  be  measured 
accurately.  With  the  measured  fringe  order  and  thickness  the  fringe 
value  2f  of  the  material  was  found  to  be  3-0  psi  per  fringe,  per  inch 
compression . 


Re  suit  s 

Interior  stresses.-  The  stress  distribution  was  obtained  along  six 


lines  passing  through  the 
cated  in  figure  21  by  the 
the  six  lines  are  defined 

A -A 

intersection  of 

B-B 

intersection  of 

C-C 

intersection  of 
to  meridian 

D-D 

load  line 

E-E 

intersection  of 
meridian 

F-F 

intersection  of 

meridian  plane  perpendicular  to  the  above  two  slices 


In  this  problem  advantage  was  taken  of  the  symmetry  of  the  sphere  to 
eliminate  the  need  for  two  models  as  discussed  in  the  theoretical  part 
of  the  report.  In  particular  because  of  the  rotational  symmetry  of  the 
stresses  one  meridian  slice  could  be  taken  to  represent  all  meridian 
slices . 

Typical  calculation.-  In  order  to  make  clear  the  application  of  the 
method  the  complete  calculations  for  line  C-C  will  now  be  given.  The 
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evaluation  of  the  stresses  along  this  line  requires  all  the  generality 
which  would  he  encountered  in  a  body  devoid  of  symmetry.  The  basic 
data  for  the  determination  of  the  stresses  on  this  line  are  obtained 
from  the  stress  patterns  and  isoclinics  at  normal  incidence  of  the  two 
slices  defining  the  line  C-C  and  from  the  stress  pattern  and  isoclinics 
of  one  of  the  slices  at  oblique  incidence  with  rotation  about  an  axis 
perpendicular  to  C-C  and  lying  in  the  plane  of  the  plate. 

Because  of  the  symmetry  of  the  stresses  along  C-C  it  is  necessary 
only  to  deal  with  half  the  length  of  the  line.  This  half  length  was 
divided  into  10  equal  subdivisions.  The  two  necessary  auxiliary  lines 
were  drawn  parallel  to  it  in  each  of  the  two  orthogonal  planes  and 
spaced  the  length  of  one  subdivision  apart. 

The  first  step  is  to  obtain  the  distribution  of  the  shearing 
stresses  along  line  C-C  and  along  the  four  auxiliary  lines.  This 
requires  the  determination  of  the  distribution  of  the  secondary  principal 
stress  differences  and  of  the  isoclinics  along  these  lines.  Figure  22 
shows  the  curves  of  n"'  and  (f for  the  slice  parallel  to  the  meridian, 
and  figure  23  shows  n"  and  0"  for  the  slice  parallel  to  the  equator. 
With  these  data  the  required  shearing  stresses  are  computed  at  each 
division  point  of  C-C.  Thus,  following  equations  (12) 

n  „  =  —  n"’  sin  20"’  ( l6a) 

y^  2 


"xz  -  \  n"  sin  ^ " 


(16b) 


It  will  be  noted  that  for  positive  values  of  z  the  shear  systems  n yZ 
and  nXz  are  positive. 


As  noted  in  the  theoretical  part  of  the  report  the  integration 
requires  the  use  of  the  difference  between  the  shearing  stresses  at  the 
center  of  each  subdivision.  These  differences  are  obtained  from  the 
curves  representing  the  distribution  of  the  shearing  stresses  just 
found.  Figure '2h  shows  the  curves  of  the  shear  differences  for  the  two 
slices . 

The  next  step  is  to  obtain  the  value  of  the  normal  stress  nz  at 

each  division  point  by  an  integration  of  one  of  the  differential  equations 
of  equilibrium.  The  integration  is  carried  out  by  approximation  using 
summations  to  replace  the  integrals.  The  appropriate  equation  for 
line  C-C  is  similar  to  equation  (lla),  that  is. 
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Choosing  Ax  =  Ay  =  Az  the  ratios  offthese  quantities  are  unity.  The 
signs  of  the  ratios  depend  on  the  choice  of  axes  and  the  direction  of 
integration.  The  integration  begins  at  the  boundary  and  proceeds  inward 
Consequently  Az  is  negative.  The  shear  differences  have  been  formed 
in  such  a  way  that  Ay  is  positive  and  Ax  negative.  The  final  form 
of  the  equation  is  therefore 

.  d  i 

(nz)j  =  Wc  +  IT  ^yz  -  £  ^xz  (17b) 


The  signs  of  the  shear  differences  are  as  shown  by  the  curves. 

In  order  to  start  the  integration  the  value  of  ^nzj  is  required. 

This  must  be  determined  from  the  boundary  conditions  and  boundary  fringe 
orders.  From  the  fact  that  the  boundary  is  unloaded  it  is  evident  that 
the  principal  stress  normal  to  the  boundary  is  zero.  Inspection  of  the 
meridian  slice  shows  that,  excepting  the  loaded  regions,  the  boundary 
stresses  in  the  meridian  section  are  also  zero.  The  fringe  order  at  the 
boundary  of  the  slice  parallel  to  the  equator  is  O.58  fringe  tension, 
and  the  direction  of  this  stress  makes  an  angle  of  55-6°  with  line  C-C. 
The  boundary  value  of  nz  is  thus  found  from  the  equations  of  stress 
transformation  as  follows: 


(nz]c  =0.58  cos255.6  =  0.185 
The  expression  for  ^nz^  therefore  takes  the  form 

j  5 

(nz)j  =  °-l85  +  55  Anyz  -  52  Anxz 


(18) 


(17c) 


The  integration  is  easily  carried  out  in  tabular  form  as  shown  in 
table  I. 

Once  nz  has  been  found  the  values  of  nx  and  n y  are  found  from 
expressions  similar  to  equation  (l3c),  that  is, 
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nx  =  nz  -  n"  cos  20" 


(19a) 


ny  =  nz  -  n"'  cos  20'”  (19b) 

■  t 

This  computation  is  shown  in  table  II. 

The  last  step  is  to  determine  the  remaining  shearing  stress 
system  nXy.  This  was  obtained  in  this  case  from  oblique  incidence 

on  the  slice  parallel  to  the  meridian  with  rotation  through  about 
the  Y-axis.  Figure  25  shows  the  fringe  order  and  isoclinic  .distribu¬ 
tions  along  C-C  for  this  case.  With  these  data  and  the  known  values 
of  nyZ  previously  determined  the  required  shearing  stress  component 

is  found  from  an  expression  similar  to  equation  (15c).  Thus 


nxy  =  \  n0y  sin  200y  -  nyz  (20) 

The  results  of  this  computation  are  shown  in  table  III.  This  completes 
the  solution  for  all  six  stress  components  along  line  C-C.  In  order 
to  determine  the  stress  components  in  pounds  per  square  inch  it  is 
necessary  only  to  multiply  the  stresses  in  fringes  by  the  proper  fringe 
value  of  the  slice. 

Using  methods  similar  to  those  just  explained  the  stress  components 
for  all  six  lines  have  been  determined.  With  the  exception  of  line  D-D 
integration  began  at  the  boundary  and  proceeded  inward.  For  line  D-D 
the  starting  point  was  taken  as  the  center  of  the  sphere  and  integration 
proceeded  upward.  The  starting  value  of  ny  for  this  line  was  taken 

to  be  that  obtained  from  line  A-A.  The  results  of  these  computations 
are  shown  in  figures  26  to  32.  At  the  center  of  the  sphere  the  stress 
components  were  found  to  be  Uy  =  -2.59P0A  and  crx  =  az  =  0.45P/A. 

These  values  may  be  compared  with  the  stresses  at  the  center  of  a  disk 
under  diametral  compression  which  are  Uy  =  -1.91P/A  and  cjx  =  0.6k'P/A. 

Contact  stresses  at  load  points.-  As  noted  previously  the  loads  on 
the  sphere  produced  considerable  local  deformation  which  resulted  in 
flattened  areas  at  the  poles.  Upon  observation  of  the  meridian  slice 
it  was  found  that  fringes  and  isoclinics  were  unusually  clear  right  to 
the  loaded  boundary.  It  was  therefore  possible  to  continue  the  integra¬ 
tion  all  the  way  to  the  loaded  boundary  along  lines  normal  to  the  areas 
of  contact  and  thereby  to  obtain  approximations  to  the  contact  stresses. 
The  stress  components  at  the  ends  of  line  D-D  represent  the  contact 
stresses  at  the  poles,  that  is,  at  the  centers  of  the  loaded  areas.  In 
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order  to  determine,  at  least  roughly,  the  actual  distribution  of  the 
contact  stresses  over  the  loaded  areas  two  additional  slices  normal  to 
the  loaded  zones  were  cut  from  the  remaining  material  of  the  sphere. 

The  intersections  of  these  slices  with  the  meridian  slice  then  define 
two  lines  parallel  to  line  D-D  which  extend  from  line  B-B  to  the  loaded 
boundary.  Starting  values  were  taken  from  the  stress  distribution  on 
line  B-B  previously  obtained,  and  integration'  proceeded  to  the  loaded 
boundary  in  the  usual  fashion.  In  this  way  the  contact  stresses  at  two 
points  at  different  distances  from  the  center  were  obtained.  With  these 
three  points  the  distribution  of  the  normal  stresses  on  the  contact 
surface  could  be  pretty  well  determined.  The  shearing  stresses  acting 
on  the  surface  of  contact  were  found  from  the  values  of  n 1  and  the 
isoclinics  0’  in  a  meridian  section  in  the  region  of  contact.  The 
results  are  shown  in  figure  33 ,  the  directions  of  the  shearing  stresses 
being  from  the  poles  outward. 

Checks  on  accuracy.-  Two  types  of  checks  are  available  in  this 
problem,  static  checks  and  checks  between  stresses  on  different  lines. 
Static  checks  were  made  from  the  stresses  on  lines  A-A  and  B-B  and  from 
those  acting  on  a  diameter  in  the  surfaces  of  contact.  Since  these 
stresses  are  rotationally  symmetrical  the  resultant  force  acting  on  the 
equatorial  plane  and  on  the  plane  containing  line  B-B  parallel  to  the 
equator  as  well  as  on  the  plane  of  contact  can  be  determined  by  integra¬ 
tion.  From  the  stresses  on  line  A-A  the  resultant  load  on  the  equatorial 
plane  was  computed  as  Tj6  pounds,  which  is  2.3  percent  higher  than  the 
applied  load  of  172  pounds.  The  stresses  on  line  B-B  gave  a  resultant 
of  168  pounds  which  is  2.3  percent  low.  Lastly  the  resultant  of  the 
normal  stresses  on  the  surface  of  contact  was  found  to  be  170  pounds, 
or  1.2  percent  low. 

Points  O’  and  Op1  (fig.  2l)  are  common  points  on  different  lines. 

The  stresses  at  these  points  can  be  determined  from  each  line  and  the 
results  compared.  The  stresses  at  O'  were  found  by  integrating  along 
line  B-B  and  also  by  integration  along  the  path  A00'.  From  line  B-B 
the  stresses  were  found  to  be  n^  =  -4.08  and  nx  =  nz  =  0.40.'  From 

the  path  A00'  they  were  found  to  be  n-y  =  -4.06  and  nx  =  nz  =  0.42. 

At  Op'  the  stresses  are  found  from  lines  C-C  and  F.-F.  From  l-ine  C-C 
the  stresses  were  computed  as  n^.  =  -1.12,  nx  =  -O.36,  nz  =  0.43,  and 

nXy  =  0.62.  From  line  F-F  they  were  n-y  =  -1.14,  nx  =  -0.39,  nz  =  0.40, 
and  nXy  =  0.66.  The  agreement  in  these  values  is  seen  to  be  quite  good. 

It  is  also  possible  to  compute  the  stresses  on  lines  C-C  and  E-E 
from  the  stresses  on  lines  A-A  and  B-B.  Figure  34  shows  the  rectangular 
stress  components  on  lines  C-C  and  B-B.  The  necessary  equations  for 
transformation  are  similar  to  the  familiar  equations  for  inclined  planes 
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in  plane  stress  systems,  that  is, 


ffy  +  CTv  (Tv  —  ffv  ,  s 

Oq  =  - - — -  +  - - — -  cos  2 0  +  Txy  sin  2 0  (21a) 


Tg  =  — g — sin  20  -  TXy  cos  20  (21b) 


The  stresses  computed  by  transformation  were  compared  with  the  stresses 
independently  determined  on  lines  C-C  and  E-E  by  integration.  The  com¬ 
parative  values  are  shown  in  figures  28,  29,  and  31*  The  general  agree¬ 
ment  is  seen  to  be  excellent. 

Check  by  the  Lame-Maxwell  equations.-  Lines  A -A  and  D-D  are  lines 
of  symmetry  for  the  sphere.  For  these  special  lines  the  stresses  can 
be  computed  by  the  method  outlined  by  Jessop  (reference  IT)  using  the 
extension  of  the  Lame-Maxwell  equations  to  three  dimensions.  This 
computation  has  been  carried  out  using  fringe  orders  from  the  meridian 
and  equatorial  slices  and  the  5°  isoclinic  from  the  meridian  slice. 

The  comparative  values  are  given  in  figure  26  for  line  A-A  and  figure  30 
for  line  D-D.  Inspection  of  the  figures  shows  that  in  both  cases  the 
agreement  is  close. 


Analysis  and  Discussion 

The  primary  objective  of  the  project  under  discussion  was  to  develop 
a  general  method  for  solving  three-dimensional  problems  photoelastic ally. 
In  the  theoretical  part  of  this  report  such  a  method  is  described.  The 
experimental  work  shows  that  the  proposed  method  is  practical. 

It  is  too  early  to  draw  broad  conclusions  regarding  the  general 
accuracy  of  the  new  method.  However,  the  excellence  of  the  static  checks 
and  the  consistency  of  the  results,  as  shown  by  the  close  cross  checks 
between  the  results  from  the  various  lines,  seem  to  indicate  possibil¬ 
ities  of  high  accuracy.  Unfortunately  there  is  as  yet  no  theoretical 
solution  available  for  this  particular  problem  to  furnish  conclusive 
checks  and  a  measure  of  the  errors.  Nevertheless  there  is  a  reasonable 
degree  of  certainty  that  the  major  stresses  are  free  from  significant 
error. 

It  must  be  pointed  out  that  the  stresses  as  found  here  represent 
the  solution  for  a  material  for  which  Poisson's  ratio  is  0.48,  whereas 
most  structural  materials  have  Poisson's  ratios  of  about  0.3-  This  is 
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an  inherent  limitation  of  three-dimensional  photoelasticity.  However, 
as  noted  previously,  the  theoretical  solutions  available  to  date  indicate 
that  Poisson's  ratio  has  only  a  small  influence  on  the  major  stresses 
although  the  effect  on  the  minor  stresses  may  be  pronounced  (refer¬ 
ences  19  and  20). 

Although  no  theoretical  solution  is  available  for  the  sphere 
Hertz's  solution  can  be  used  to  check  the  contact  stresses  determined 
photoelastically.  According  to  Hertz's  solution  the  contact  pressure 
is  representable  by  the  ordinates  to  a  hemisphere  erected  on  the  contact 
surface.  Further  the  Hertz  theory  predicts  that  the  maximum  contact 
pressure  should  be  1.5P/AC.  Reference  to  figure  33  shows  that  the 
experimentally  determined  value  of  the  maximum  pressure  is  1.53P/Ac, 
which  is  2  percent  high.  When  it  is  considered  that  the  path  of  integra¬ 
tion  used  in  determining  the  maximum  pressure  led  along  an  equatorial 
radius  to  the  center  and  thence  up  the  load  axis  to  the  surface  this  is 
indeed  a  remarkable  check.  The  general  distribution  of  the  pressures 
is  also  seen  to  be  substantially  correct.  The  proposed  method  would 
therefore  seem  to  hold  considerable  promise  for  the  determination  of 
contact  stresses. 


According  to  the  Hertz  theory  the  two  principal  stresses  in  the 
plane  of  the  contact  surface  at  the  pole,  ax  and  crz,  should  each 

equal  — --v  a  .  For  Poisson's  ratio  of  l/2  this  reduces  to 
2  * 

ax  =  az  ~  °y  which  indicates  that  an  isotropic  point  exists  at  the  pole, 

and  therefore  the  shearing  stress  is  zero  at  this  point.  This  is  borne 
out  by  the  photoelastic  results. 


Along  the  circular  boundary  of  the  surface  of  contact  the  stress 

consists  of  a  pure  shear  of  the  amount 

Hertz  theory.  If  Poisson's  ratio  is  l/2  this  expression  vanishes. 
Hence  the  normal  and  shearing  stresses  at  the  boundary  should  be  zero. 
The  photoelastic  results  are  seen  to  give  this  value. 


according  to  the 


The  problem  treated  in  this  report  has  complete  rotational  symmetry 
which  simplifies  the  experimental  technique  by  eliminating  the  use  of 
two  models.  No  problem  has  as  yet  been  solved  which  requires  two  models. 
The  use  of  two  models  will  no  doubt  introduce  complications,  but  no 
insurmountable  difficulties  are  anticipated.  However,  further  work  must 
be  done  to  demonstrate  the  effectiveness  of  the  proposed  method  with  two 
models . 


There  remains  also  the  possibility  of  using  sub  slices  from  the 
main  slice  as  discussed  in  the  theoretical  part.  The  model  used  in  this 
investigation  was  not  large  enough  to  make  this  procedure  feasible 
although  some  attempts  were  made.  This  possibility  also  needs  to  be 
further  explored. 
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It  will  perhaps  also  he  desirable  to  repeat  the  solution  of  the 
sphere  with  smaller  loads  in  order  to  reduce  the  relatively  large  local 
deformations . 


SUMMARY  OF  RESULTS 


The  results  from  this  investigation  to  develop  a  general  method  for 
three-dimensional  photoelastic  stress  analysis  may  be  summarized  as 
follows: 

1.  The  method  of  strain  measurement  after  annealing  cannot  be  used 
with  the  materials  now  available  in  this  country. 

2.  A  general  photoelastic  method  of  obtaining  six  stress  components 
at  any  point  of  an  unsymmetrical  body  arbitrarily  loaded  has  been 
developed.  This  method  does  not  depend  on  Poisson's  ratio,  although  the 
results  reflect  the  physical  constants  of  the  model. 

3.  The  new  method  is  applicable  in  the  plastic  range  of  the  model. 

4.  The  method  shows  possibilities  for  the  determination  of  contact 
stresses . 

5.  The  stresses  existing  in  a  sphere  subjected  to  diametral  com¬ 
pression  have  been  determined  with  considerable  accuracy. 

6.  At  the  center  of  the  sphere  the  stress  components  were  found  to 
be  cry  =  -2.59P/A  and  ax  =  oz  =  0.45P/A,  where  P  is  the  load  on  the 

sphere  and  A  is  the  area  of  the  equatorial  plane  of  the  sphere.  These 
values  may  be  compared  with  the  stresses  at  the  center  of  a  disk  under 
diametral  compression  which  are  Oy  =  -1.91P/A  and  crx  =  0.64P/A. 

7.  Further  work  is  needed  to  determine  the  full  potentialities  of 
the  method  when  two  models  are  used.  Further  work  is  also  needed  to 
determine  the  practicability  of  sub  slices. 


Illinois  Institue  of  Technology 

Chicago  16,  Ill.,  September  15,  1951 
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TABLE  III 

CALCULATION  OF  SHEARING  STRESS  BY  OBLIQUE  INCIDENCE  FOR  LINE  C-C 
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Figure  2.-  Necessary  slices  and  directions  of  light  for  determination  of 
surface  stresses.  L,  direction  of  light. 
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Figure  8.-  Stress  pattern' of  whole  sphere. 
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Figure  9*-  Slicing  plan  for  sphere.  1,  slice  parallel  to  equator; 
2,  meridian  slice;  3 >  slice  parallel  to  meridian;  and  k,  equa¬ 
torial  slice. 


Figure  15.-  Stress  pattern  of  slice  parallel  to  meridian  for  a  rotation 

of  45°  about  Y-axis. 
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Figure  l6.-  Typical  isoclinic  for  meridian  slice. 


Figure  17. -  Typical  isoclinic  for  slice  parallel  to  meridian 


Figure  19.-  Stress  pattern  of  calibration  member. 
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Figure  21.-  Lines  in  sphere  along  which  stress  distributions  were 

determined. 
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eter  <p' '  '  for  line  C-C  and  wo  auxiliary  lines  m  plice  para, 
to  meridian.  In  curves  II  and  III  the  letter  C  denotes  the 
X-coordinate  of  line  C-C. 
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Figure  28.-  Distribution  of  normal  stresses  along  line  C-C. 


Dimensionless  Stress, 
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Figure  30.-  Distribution  of  stresses  along  line  D-D. 
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